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Abstract 



We consider the production of a generic system of non-strongly interacting particles 
with a high total invariant mass M in hadron collisions. We examine the transverse- 
momentum (gy) distribution of the system in the small-gr region (q? <C M), and we 
present a study of the perturbative QCD contributions that are enhanced by powers 
of large logarithmic terms of the type ln(M 2 /gf.). These terms can be resummed 
to all orders in QCD perturbation theory. The partonic production mechanism of 
the final-state system can be controlled by quark-antiquark (qq) annihilation and/or 
by gluon fusion. The resummation formalism for the qq annihilation subprocess is 
well established, and it is usually extrapolated to the gluon fusion subprocess. We 
point out that this naive extrapolation is not correct, and we present the all-order 
resummation formula for the qr distribution in gluon fusion processes. The gluon 
fusion resummation formula has a richer structure than the resummation formula in 
qq annihilation. The additional structure originates from collinear correlations that are 
a specific feature of the evolution of the colliding hadrons into gluon partonic states. 
In the qT cross section at small values of qx, these gluon collinear correlations produce 
coherent spin correlations between the helicity states of the initial-state gluons and 
definite azimuthal-angle correlations between the final-state particles of the observed 
high-mass system. 
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1 Introduction 



The properties of the transverse-momentum distributions of systems of high invariant mass that 
are produced at high-energy hadron colliders are important for QCD and electroweak studies and 
for physics studies beyond the Standard Model (SM). 

We consider the inclusive hard-scattering process 

h 1 (p 1 ) + h 2 (p 2 )^F{{q i })+X , (1) 

where the collision of the two hadrons hi and h 2 with momenta p\ and p 2 produces the triggered 
final-state system F, accompanied by an arbitrary and undetected final state X. We denote by 
sfs the centre-of-mass energy of the colliding hadrons, which are treated as massless particles 
(s = (pi + p 2 ) 2 = 2pip 2 ). The observed final state F is a generic system of one or more particles 
with momenta gf (i = 3, 4, 5, . . . ). The total momentum of F is denoted by q^ (q = J2i 9i)> an d it 
can be expressed in terms of the total invariant mass M (q 2 = M 2 ), the transverse momentum qx 
with respect to the direction of the colliding hadrons, and the rapidity y (2y = \n(p 2 q/piq)) in the 
centre-of-mass system of the collision. Throughout the paper, we limit ourselves to considering 
the case in which the system F is formed by non-strongly interacting particles, such as vector 
bosons (7, W, Z, . . . ), Drell-Yan (DY) lepton pairs, Higgs particles and so forth. 

Provided the invariant mass M is large (M 3> Aqcd, ^qcd being the QCD scale), the 
production cross section and associated kinematical distributions of the process in Eq. (CQ), can be 
evaluated by using QCD perturbation theory. The cross section is expressed as a convolution of 
partonic cross sections da a b (a, b = q,q, g) with the parton densities of the colliding hadrons. The 
partonic cross section da a b is computed as a power series expansion in the QCD coupling a$(M 2 ) 
by considering the corresponding partonic subprocess a + b —> F + . . . , where the dots denote 
final-state partons. Since F is a system of colourless particles, the partonic subprocesses include 
quark- ant iquark (qq) annihilation, 

q + q -> F , (2) 

and gluon fusion, 

g + g^F . (3) 

In this paper we are interested in considering the process of Eq. (CQ) in kinematical configu- 
rations where the transverse momentum qx of the system F is small (say, gr M). Unless the 
subprocesses in Eqs. (T5]) and ([3]) are both forbidden by selection rules related to the nature of F 
(e.g., gg — >■ W ± is forbidden), these subprocesses produce the system F with q T = 0. The system 
F acquires a non- vanishing transverse momentum through higher-order QCD radiative corrections 
to the subprocesses in Eqs. (j2j) and ([3]). Nonetheless, the bulk of the events is still produced in the 
small-gr region. 

The perturbative-QCD computation of the partonic cross sections da a b in powers of as(M 2 ) 
shows that high-order coefficients contain logarithmic terms of the type ln m (M 2 /g|,). Although 
as(M 2 ) is small, these logarithmic terms can be large in the small-gr region (qx <C M), thus 
spoiling the quantitative convergence of the expansion in powers of as(M 2 ) (at each fixed order in 
as, the partonic cross section eventually diverges to either +00 or —00 by considering the limit 



1 



qr — > 0). To obtain reliable perturbative predictions in the small-g-r region, the logarithmically- 
enhanced terms have to be evaluated at sufficiently-high perturbative order^j], and possibly re- 
summed to all orders in as(M 2 ). 

The small-g^ logarithmic terms have their physical origin from multiple radiation of final-state 
partons that are soft and/or collinear to the colliding hadrons (partons). The method and the 
formalism to resum the logarithmically-enhanced terms at small qr was developed in the eighties 
[TT-fTT], Subsequent, and important, theoretical progress in this field regards, for instance, the 
explicit computation of high-order resummation coefficients [T2HT5] (see additional comments in 
Sect. [2]) and the understanding of their universality (process-independent) structure [T2| [T6]. 

In this paper we consider small-gr resummation, and we deal with an issue that is also related to 
universality. The issue regards the relation between processes that are controlled by qq annihilation 
and by gluon fusion. Transverse-momentum resummation was originally worked out for the DY 
process fTHTO] , which is driven by qq annihilation. The resummation structure that emerges in the 
DY process was then customarily used (see, e.g., Refs. [T7H26] and references therein) for many 
other processes of the class in Eq. ([T]). Such processes include, for instance, the production of the 
SM Higgs boson [TT1 - I2T] . of photon pairs [22], of vector boson pairs such as ZZ [23] and W + W~ 
[23], and of slepton pairs [251 EE]- 111 particular, Higgs boson production is driven by gluon fusion, 
whereas diphoton and diboson production receives contributions from both the qq annihilation 
and gluon fusion subprocesses. 

In the present contribution, we point out that there are some key differences between qq 
annihilation and gluon fusion. These differences have escaped detection until recent findings [T3| 
1271 128] , The physical origin of the differences is due to specific collinear correlations (see Sect. [3]) 
that are a distinctive feature of the perturbative evolution of the colliding hadrons into gluon initial 
states. Analogous correlations are not produced by the perturbative evolution of spin unpolarized 
hadrons into quark or antiquark initial states. 

As a consequence of these differences, transverse-momentum resummation in gluon fusion sub- 
processes has a 'richer' structure than in qq annihilation subprocesses. The small-g^ resummation 
formalism for the DY process [THTU] has to be modified and extended to deal with gluon fusion 
subprocesses. In particular, in gluon fusion subprocesses, gluon collinear correlations produce spin 
and azimuthal correlations that are logarithmically enhanced in the small-gr region. 

In the following we present and discuss our main general results on transverse-momentum 
resummation in gluon fusion processes. Then, we shall comment on Refs. [131 EZl EE]. Details 
about the derivation of our results, and the illustration of further related results, will appear 
elsewhere [2§] . 

The outline of the paper is as follows. In Sect. |2J we briefly review the classical QCD results 
on transverse-momentum resummation. These results, which mostly derive from studies of the 
DY process, are presented using a general and process-independent notation that is useful for the 
subsequent presentation of qx resummation in gluon fusion processes. In Sect. [31 we present our 
all-order resummation formula for generic transverse-momentum cross sections controlled by gluon 
fusion. We illustrate the structure of the resummation formula, and we discuss its origin from 

^The 'sufficiently-high' order depends on the specific qx region of interest in each specific process; this order 
cannot be specified 'a priori'. 
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quantum-mechanical correlations (interference effects) produced by the collinear-parton radiation 
that accompanies the gluon fusion hard-scattering subprocess. We also explicitly consider the 
specific example of SM Higgs boson production. In Sect. HI we reformulate qx resummation in 
the helicity space of the colliding gluons. We show how gluon collinear correlations are related to 
helicity-flip phenomena in the hard-scattering subprocess. In Sect. El we specify the gluon fusion 
resummation formula for azimuthally-averaged transverse-momentum cross sections. We point 
out that the differences between the qq annihilation and gluon fusion channels persist even after 
having performed the integration over the azimuthal angle of the transverse-momentum vector. 
Section |6] is devoted to derive and discuss the general structure of the azimuthal-angle correlations 
embodied in the gluon fusion resummation formula. Few summarizing remarks are presented in 
Sect, m 



2 Small-g^ resummation in impact parameter space 



In this section we recall the 'classical' formalism [THTT] [16] of transverse-momentum resummation 
in impact parameter space. This illustration sets the stage for the presentation of our results on 
small-gr resummation in gluon fusion processes (see Sects. EHI!])- 

We consider the process in Eq. (0Q), and we introduce the corresponding multidifferential cross 
section 

d^d^dydCl (P^T,M,y,tl) . (4) 

The differential cross section depends on the total momentum of the system F (i.e. on the variables 
qx,M, y) and, to be quite general, it can also depend on additional variables that specify the 
kinematics of the particles in the system F. In Eq. (j4|) these additional variables are generically 
denoted as ft = {Qa, fie, • • • } (correspondingly, we define dfl = cffi^dOs . . . ). They can be, for 
instance, the rapidity yi and the azimuthal angle 0(qxj of one of the particles (with momentum 
qi) in the system F . In general, we only assume that the kinematical variables {f^, fig, ... } are 
independent of qx, M and y. 

Considering the qx dependence of the multidifferential cross section in Eq. (j3J) within pertur- 
bative QCD, we introduce the following decomposition: 

da F = daf ng) + dap s) . (5) 

Both terms on the right-hand side are obtained as convolutions of partonic cross sections and the 
scale-dependent parton distributions f a /h{x,n 2 ) (a = q/,q~f,g is the parton label) of the colliding 
hadronq3 The distinction between the two terms is purely theoretical. The partonic cross sec- 
tions that enter the singular component (the first term on the right-hand side) contain all the 
contributions that are enhanced (or 'singular') at small q?. These contributions are proportional 
to <^ 2 )(qx) or to large logarithms^ of the type \ ln m (M 2 /g|,). On the contrary, the partonic cross 



* Throughout the paper we always use parton densities as defined in the MS factorization scheme, and as(q 2 ) is 
the QCD running coupling in the MS renormalization scheme. 

§To be precise, the logarithms are combined with corresponding 'contact' terms, which are proportional to 
<5( 2 )(qx). These combinations define regularized (integrable) 'plus distributions' t- ln m (Af 2 /g£) 



to q T . 



with respect 

+ 
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sections of the second term on the right-hand side are regular (i.e. free of logarithmic terms) 
order- by-order in perturbation theory as qx — > 0. To be precise, the integration of dap Cg ^ /d 2 q^r 
over the range < qx < Qo leads to a finite result that, at each fixed order in a$, vanishes in the 
limit Qo -> 0. 

The regular component da F eg ^ of the qx cross section is definitely process dependent. In this 
paper we limit ourselves to considering the singular component, which has a universal (process- 
independent) structure. 

To simplify the presentation of the all-order (resummed) structure of the singular component 
of the qx differential cross section in Eq. (j3j), we introduce a shorthand (symbolical) notation in 
several places. For instance, the singular component of Eq. (TJJ is simply denoted by [dap]; namely, 
we define 

Ax (sing) 

[dap] S <Pq T dM 2 dy dn {PlM qT ' M ' n) ■ (6) 
The transverse-momentum resummation formula can be written in the following factorized form 

mm. 



c=q,q,g 

£ f— f— [^^] c , aia2 f ai / hl fa/zi, b 2 /b 2 ) f a2/h2 (x 2 /z 2 , bl/b 2 ) , (7) 



X 

ai,a2 



where bo = 2e lE (7^ = 0.5772 ... is the Euler number) is a numerical coefficient, and the 
kinematical variables x\ and x 2 are 

X! = ^= e +y , x 2 = ^= e~y . (8) 

V-S y/S 

The function S C (M, b) and the functions symbolically denoted by da'P and [H F CiC 2 \ are 
specified below. 

The right-hand side of Eq. ([7]) involves the Fourier transformation with respect to the impact 
parameter b and two convolutions over the longitudinal-momentum fractions Z\ and z 2 . The 
parton densities f a ./ h .(x, fi 2 ) of the colliding hadrons are evaluated at the scale \i = bo/b, which 
depends on the impact parameter. 

We note that in the context of the study of the present paper, the resummation formula in 
Eq. (JTj), and the resummation formulae in Sects. [3HSI have a purely perturbative-QCD contend, 
analogously to customary fixed-order calculations of hard-scattering cross sections in hadron col- 
lisions. Using the Altarelli-Parisi evolution equations, the parton densities f a /h{x,bl/b 2 ) can be 
expressed [TJ] in terms of the corresponding parton densities f a /h(x,fi 2 F ) at the evolution scale 
\i — fip, where /xp is the customary factorization scale that enters fixed-order calculations. Hav- 
ing done that, all the remaining factors on the right-hand side of Eq. (j7]) are partonic contributions 
that can be expanded in powers of as(M 2 ) at arbitrary perturbative orders. 



^Throughout the paper we do not consider the inclusion of any non-perturbative contributions, such as, for 
instance, those first introduced in Rcf. |5j. 
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The small-g-r region where C M corresponds in impact parameter space to the large-6 region 
where b 3> 1/M. The perturbative expansion of the b space integrand in Eq. (j7]) produces large 
perturbative coefficients of the type ln m (6 2 M 2 ): these coefficients lead to the small-g^ logarithmic 



terms 



4 In" 1-1 
q T space. 



(M 2 /q 



through the evaluation of the Fourier transformation from b space to 



The factor 



da 



(o) 



cc,F 



in Eq. (J7|) depends on the process (i.e. on the specific final state system 

F and its kinematics). This factor is the Born level cross section da^ (i.e. the cross section at 
its corresponding lowest order in as) of the partonic subprocesses c + c — > F in Eqs. (J2J) and (j3J). 
Making the symbolic notation explicit, we write: 



da 



(0) 
cc,F 



da 



(0) 
cc,F 



M 2 dfl 



(9) 



where X\p^ {X2P2) is the momentum of the parton c (c). In Eq. (J7J), we have included the contri- 
bution of both the qq annihilation^] channel (c = q, q) and the gluon fusion channel (c = g); one 



of these two contributing channels may be absent (i.e. 
on the specific final state F. 



da 



(0) 
cc, F 



in that channel), depending 



The factor S c (M,b) in Eq. (J7|) is universal (process independent): it does not depend on the 
produced final-state system F and on its kinematics. It only depends on the partonic channel that 



produces the cross section 



da 



(0) 
cc, F 



Sudakov form factor in the cases o 



. Thus, S C (M, b) is called quark (c = q or q) or gluon (c = g) 
the qq annihilation or gluon fusion channel, respectively. The 



Sudakov form factor can be expressed in the following exponential form: 



M 



A c (a s (q 2 )) \n— + B c (a s (q 2 )) 



S c (M,b)=exp\- [ % 

[ Jbyv q 2 

where the functions A c (as) and B c (as) are perturbative series in a s : 

00 

Ma.) - £(?)>>• 

71=1 

OO 



(10) 



n=l 



ID 



(12) 



The factor \H F CiC2 1 \ in Eq. (JTJ) has the following explicit form: 

[H F C 1 C 2 ] c5ma2 =H F (x lPl ,x 2 p 2 ;n;a s (M 2 )) C cai ( Zl ^s(bl/b 2 )) C 5a2 (z 2 ;a s (b 2 /b 2 )) , (13) 
where H F and C a b are both functions of as, and they can be perturbatively expanded as follows: 



00 

H[(x lPl ,X2P2;^a s ) = 1 + J2(^Y H[ {n) (x lPl ,X2P2;n) 

n=l 

00 

C ab (z;a s ) = 5 ab 5(l-z) + J2(^yC i : b ) (z). 



(14) 
(15) 



71=1 



"in the case of qq annihilation, the notation c = q, q is not completely precise, since the quark and the antiquark 
can have either equal or different flavour. The same comment applies to the factor [.ff^CiC^] cg . Q a ■ 
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The function is process dependent, whereas the perturbative functions C a b are universal 
and depend only on the parton indices a and b (analogously to the dependence of the anoma- 
lous dimensions 7 a b that control the perturbative evolution of the parton densities, through the 
Altarelli-Parisi equations). 

By inspection of the right-hand side of Eq. ( IT3l) . we notice that the scale of as is not set to 
a unique value. We have as(M 2 ) in the case of the function iff (as naturally expected for a 
process-dependent contribution), and as(&o/& 2 ) in the case of the functions C cai and C 5a2 . The 
replacement as(M 2 ) — > as(&o/& 2 ) in Eq. f|T3|) is feasible, provided it is properly compensated [16] 
by a corresponding factor to be inserted in Eq. ffTUl) : this procedure leads to a modification of the 
Sudakov form factor, which becomes a process-dependent quantity. Indeed, the hard-scattering 
function was introduced in Ref. [IS] to explicitly show (see also Ref. [12]) the universality 
(process independence) of both the Sudakov form factor and the coefficient function C a b. In the 
version of Eq. ([7]) that was originally presented for the DY process (see Eq. (1.1) in Ref. [TO]), 
the function = H® Y (note that H® Y depends only on as) is absorbed in the definition of the 
functions C a b and of the function B c of the form factor (these functions are thus 'those' of the 
DY process). As shown in Sect. EJ the process-dependent hard-scattering function definitely 
plays a distinctive role in the case of gluon fusion subprocesses. 

The present knowledge of the perturbative coefficients in Eqs. (TTTj) . (jT2j) . (|T^|) and (|T5l) is as 
follows (see also Sect. 2.3 of Ref. [20] . where we used the same notation as in the present paper). 
The coefficients A { c l \ and A c 2) are known since a long time for both the quark [6] and the 
gluon p~T] form factors. The explicit expression of the coefficient for the DY process was first 
presented in Ref. [9]. The process- independent structure and the explicit form of the coefficients 
Be (c = q, g) was derived in Ref. [12] . The result of A c has been obtained very recently p~5] 
by relating its value to the coefficient Be 0, [12] and to the coefficient of the soft part of the 
Altarelli-Parisi splitting functions at 0(a%) |30] . The universal first-order coefficients Cqg (z) and 
Cg^z) were first computed in Refs. [H [10] and Ref. [51] . respectively. The general result for the 
first-order coefficients He and C^ h }{z) was derived in Ref. [12] . where the process dependence 

F (1) 

of H c is explicitly related to the first-order virtual corrections of the partonic subprocesses 
c + c — > F in Eqs. ([2]) and ([3J). The coefficients of Eq. ([TBI at the second order in as have been 
computed for both SM Higgs boson production by gluon fusiorj**! [T3] and the DY process [Ml [32] , 

We note that the resummation formula ([7]) involves perturbative functions of as(/i 2 ), where the 
scale /i 2 is /i 2 = M 2 (see Eqs. and (QJ)), or fi 2 = q 2 (see Eq. (flOl)), or /i 2 = b 2 /b 2 (see Eq. (fT3l)). 
All these functions can be expressed in terms of as(fi 2 R ), ln(/z 2 / /z 2 ^) and the perturbative coefficients 
of the QCD /3-function (/ir is the renormalization scale that customarily appears in fixed-order 
calculations) by using the renormalization group equation for the perturbative evolution of the 
QCD running coupling as(/x 2 ). 

We add a relevant (though known) observation. Considering the dependence on the impact 
parameter b, all the factors in the integrand of the Fourier transformation on the right-hand side 
of the resummation formula ([7]) are functions of b 2 , with no dependence on the azimuthal angle 
0(b) of b in the transverse plane of the collision. Therefore, in Eq. (J7|) we can straightforwardly 

**See, however, related comments in Sect.O 
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perform the integration over 0(b) and implement the replacement 

J eJb qT ^( b2 ) = J dbb J (bq T ) F{b 2 ) , (16) 

where Jq(x) is the Oth-order Bessel function, and F(b 2 ) denotes a generic function of b 2 . This 
result implies a technical simplification of the resummation formula, since the iwo-dimensional 
Fourier transformation is replaced by the one-dimensional Bessel transformation. More impor- 
tantly, this implies that the right-hand side of Eq. ([7j) depends only on q T 2 , with no additional 
dependence on the azimuthal angle 0(q-r) of q T . Therefore, according to Eqs. (j7]) and ffTB"]) . the 
singular part of the qx differential cross section in Eq. ([6]) does not contain any azimuthal cor- 
relations with respect to qx- Equivalently, we can say that dopld 2 ^ and dopjdq 2 -. are simply 
proportional, and we can write: 

A (sing) 1 ^( sin s) 
aa F _ £ aa F 

d 2 q T dM 2 dy dCt ir dq 2 dM 2 dy dQ ' 1 ' 

Obviously, this does not mean that the multidifferential cross section dapjd 2 ^ in Eq. (jlj) has no 

4 



azimuthal correlations. In general, azimuthal correlations are present in the regular part daf eg! 



(see Eq. ([5])) of the cross section. In the small-gy region (qp -C M), dap^^ is of 0{qT/M) (modulo 
powers of ln(M 2 /q T )) with respect to d<jp mg ^ order-by-order in QCD perturbation theory. In the 
case of gluon fusion processes, these conclusions about azimuthal correlations at small qp are no 
longer true in view of the results presented in the next section. 



3 Transverse- momentum resummation in gluon fusion 

processes 

The resummation formalism reviewed in Sect. [2] was originally developed and proven [THTU] for 
the DY process (and related observables, such as the energy-energy correlation function in e + e~ 
annihilation). In the subsequent literature these results were extrapolated to various processes of 
the class in Eq. ([1]). In this section, we show that this 'naive' extrapolation is not valid in the case 
of the (sub)processes that are controlled by gluon fusion. 

To present our results, we start from the resummation formula in Eq. (J7J). It includes the 
contributions from both qq annihilation (c = q, q) and gluon fusion (c = g). We thus separate 
these two types of contributions, and we write: 

[da F ] = [da F ] {q " n - ] + [dap]^^ . (18) 

Our new results refer to [dap] ins '\ whereas we maintain the results of Sect. [2] for [dap\^ q ^ &nn '\ 
To be precise, also in the case of gluon fusion, we confirm the factorization structure on the 
right-hand side of Eq. ([7j). We explicitly report this structure: 



M 



2 



S 



d 2 b 



(2tt. 



01,02 



/ - 1 / - 1 \^ F G X C^ f ai/hl (x 1 /z 1 ,b 2 /b 2 )f a2/h2 (x2/z2,bl/b 2 ) ,(19) 
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where the lowest-order cross section 
Eqs. © and ffTUj) . respectively. The new resn 



and the gluon form factor S g (M,b) are given in 
ts regard the gluon fusion factor \H F CiC 2 ] 



The naive expression on the right-hand side of Eq. ( TT3i) has to be replaced by the following 
result (29]: 

l HFC ^] ggm = HF iuiii2U2 ( XlPl , X2P2 ;n;a s (M*)) 

x C^(z 1 ;p 1 , P 2,b;a s (b 2 /b 2 )) C^(z 2 ; Pl ,p 2 ,b;a s (b 2 /b 2 )) . (20) 

The first evident difference with respect to Eq. (1131) is the presence of Lorentz tensors, rather 
than scalar functions. The Lorentz indices (symbolically) refer to the gluon fusion hard-scattering 
process 

9(m)(xiPi) + 9{n2){x2V2) -> F , (21) 

where m (i — 1, 2) is the Lorentz index carried by the external gluon leg with incoming momentum 
XiPi. The indices v\ and v 2 refer to the external gluon legs of the process that is complex conjugate 
to that in Eq. fl2T]) . 

The pro cess- dependent factor H F in Eq. ( 120]) has the following perturbative expansion: 

oo 

+ E(v) n ^ ( ^^ (xiPi ' X2P2;f2) ' (22) 



7T 

n=l 



with the lowest-order constraint: 

H X m 9^ vl g^ = l • (23) 



We also define the scalar function H F as follows: 

H F (x lPl ,x 2 p 2 ; n, a s ) = H F % ^ V2 (x lPl , x 2 p 2 ; Q, a s ) g>* ^ g^ v * . (24) 

A relevant property [29] of the Lorentz tensor H F Ui M2 U2 is current conservation, namely: 

Pi 1 H F 11/lilX2l/2 (x lPl ,x 2 p 2 ;Sl;a s ) =p u 1 1 H F 1Vl ^ 2U2 (x lPl ,x 2 p 2 ;fl;as) = , 

P? H^ lVl ^ 2V2 {xipi,x 2 p 2 ]Q,- 1 as) = p 2 2 H F lVu/X2l/2 (x lPl , x 2 p 2 ;fl; as) = . (25) 



The universal (process-independent) partonic tensor C^ iUi in Eq. (1201) exhibits an explicit 
dependence on the impact parameter b, besides the implicit dependence on b 2 through the scale 
of as- The structure of the partonic tensor is: 

C g ^(z; Pl ,p 2 , b; as) = d^( Pl ,p 2 ) C ga (z; a s ) + D^( Ph p 2 ; b) G ga (z; a s ) , (26) 

where 

d^( P uP2) = -g^ + P ' pU2+P * pUl , (27) 

Pi -P2 

D^( Pl ,p 2] b) = d^( Pl ,p 2 ) -2 b -^- , (28) 



and b^ = (0, b, 0) is the two-dimensional impact parameter vector in the four- dimensional notation 
(b^b^ = — b 2 ). The gluonic coefficient function C ga (z;as) has the same perturbative structure 
as in Eq. ( |T5l) . The first-order coefficient C^l{z) of the function C ga (z;as) in Eq. ( [26]) and the 
first-order coefficient Hg of the function Hg in Eq. ([24"]) have actually the same value [12] as 
obtained in the context of the 'naive' expression in Eq. (|T3|) . 

The partonic coefficient function G ga (z; as) in Eq. (126]) is a specific and distinctive feature of 
transverse-momentum resummation in gluon fusion processes^. Its perturbative expansion starts 
at order as, and we write: 

oo 

G ga (z;a s ) = ^ G«(s) +E (v)"^^ • (29) 

n=2 

The first-order coefficients are [29] 

Gfl(z) = C A ^ , (30) 
G g 1 l{z) = G g 1 l{z) = C F — . (31) 

The second-order coefficients CgJ(z) and Gpo(^) are considered in Refs. [TBI 129], 

The tensors in Eqs. (127|) and ([28]) fulfil the relations g^d^ = —2 and g^ v D^ v = 0. Considering 
the centre-of-mass system of the collision and denoting by fi = 1,2 the Lorentz indices of the 
non- vanishing components of purely-transverse vectors (such as qx and b), we have = D^ v = 
if // = (v = 0) or // = 3 (v = 3). Therefore, the only non- vanishing components of d^ u and D^ v 
are those that correspond to Lorentz indices j = 1,2 and k = 1, 2 of the transverse plane; we have 

d 3k ( P i,P2) = -g 3k , D^(p uP2 ;b) = -g^-2^ . (32) 



The structure of Eqs. (120]) and (126]) has a definite physical origin: this structure is produced by 
gluon collinear correlations [29] . We refer to the correlations that occur in the universal (process- 
independent) partonic subprocess 

a ->• g + ax + a 2 + . . . , (33) 

where the initial-state colliding parton a (a = q, q, g) 'evolves' in the colliding gluon g through 
collinear radiation of the final-state partons a 1; a 2 , ■ ■ ■ ■ 

To illustrate the role of gluon collinear correlations, we briefly sketch how they arise at the first 
non-trivial order in QCD perturbation theory. We consider the partonic hard-scattering process 

a{p)+g(p) F(q) + a(k) , (34) 

which leads to the first-order (real) radiative corrections to the gluon fusion process in Eq. (J3]), 
namely, 

g(p)+g(p) -> F(g) . (35) 

"Setting G ga (z; as) — in Eq. (|26l) . the gluon fusion factor of Eq. ([20]) coincides with the corresponding 'naive' 
factor of Eq. ([131) . 
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Here the parton momenta are denoted by p,p and k. In the process (I35p . the final-state system 
F is produced with a vanishing transverse momentum. The final-state system F acquires a non- 
vanishing transverse momentum (q^qr^i = — <1t 2 ) through the radiative process in Eq. (134]) ; 
here, the final-state parton a has transverse momentum hj,, and momentum conservation implies 
fcj, = — qj,. In the small-gr region (formally, when qx — > 0), the scattering amplitude of the 
process in Eq. AMI) is singular, and the singular behaviour is controlled by a well-known universal 
factorization formula (see, e.g., Eq. (4.23) and related formulae in Sect. 4.3 of Ref. [33]). The 
factorization formula relates the scattering amplitude of the process in Eq. (|34|) with the scattering 
amplitude of the hard-scattering subprocess g(zp) + g(p) — » F. The singular factor is produced 
by the transverse-momentum spectrum of the collinear splitting subprocess 



a(p) -> g(zp) + o((l - z)p) 



(36) 



where z is the longitudinal-momentum fraction that is transferred from the initial-state parton a 
to the initial-state gluon g. At the lowest perturbative order in the QCD coupling, the transverse 
momentum spectrum of the subprocess in Eq. fl36|) is proportional to 



«s d 2 q T 

5- dz 

it q T 2 



>J5(*,qT) 



(37) 



where ktp = — qi^ is the transverse momentum of the the final-state parton a in Eq. (136]) . and the 



collinear splitting function 



has the following explicit from: 



P$(z,<tr) 



ft V 



2C A 



-g 



fj, V 



1 - Z 



+ z(l-z))+2 



1 1 n 



2C f 



q^q^ 1 — z 
q T 2 z 

,^ 1 -z + 2 q M — 



q T z 



z 



(38) 



(39) 



In Eq. ( )37|) . /x and z/ are the Lorentz indices of the gluon g(zp) in the process (136|) and in its 
complex conjugate process, respectively. 

In the context of our present study, the most important feature of the gluonic splitting process 
in Eq. (|3"6"|) is that it is intrinsically polarized. The corresponding collinear splitting function 

Pga(z, q T ) has a non-trivial dependence on the Lorentz (and, thus, spin) indices of the gluon, 
and this dependence is controlled by the azimuthal angle of the small transverse momentum that 
is radiated in the splitting process. We remark that this polarization effect is present despite 
the fact that we are not considering polarized-scattering processes. Indeed, we have performed 
the sum over the spin polarizations of the final-state parton a((l — z)p) and the average over 
the spin polarizations of the initial-state parton a(p). The intrinsic gluon polarization effects 
that arise in the splitting process (I3"6"|) (and, more generally, in the collinear splitting process of 
Eq. (I3"3"j) ) produce correlations between the initial-state gluon legs of the scattering amplitude of 
the factorized hard-scattering subprocess g(zp) + g(p) — > F and of the corresponding complex 
conjugate scattering amplitude. 

Replacing the gluon g with a quark q (or an antiquark q) in Eq. (133]) . we obtain the quark 
(antiquark) collinear evolution process 



a ->■ q(q) + a l + a 2 + 



(40) 
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In this case, if we sum over the spin polarizations of the final-state partons and we average 
over the spin polarizations of the initial-state parton a, the collinear evolution of the quark (an- 
tiquark) turns out to be unpolarized, as a consequence of helicity conservation in QCD radiation 
from a massless quark (antiquark). This essential difference between the collinear evolution of 
quarks (antiquarks) and gluons is eventually the origin of the difference of transverse-momentum 
resummation between qq annihilation processes and gluon fusion processes. 



Going back to Eqs. fl37|) - fl39|) . we can rewrite the gluon splitting functions as follows: 

Pj^qT)]" = -9"" PS{z)+^ + 2 q ^j2G^) , a = g,q,q, (41) 



where the functions Ggl(z) are those in Eqs. (130]) (a = g) and fl31~]) (a = q,q), and the functions 



Pga(z) are 



P^{z) = 2C A 



1 — z , . 
+ + z(l - z) 



1-z 

1 + (1 - z) 



(42) 



We (z) are the (real part of) first-order Altarelli-Parisi probabilities that control the 

customary collinear evolution of the unpolarized gluon parton density f g /h(%, IJ 2 )- In Eq. (j4Tp . 
the unpolarized splitting functions Pg^z) are multiplied by the tensor^ — g^ v that does not 
produce any gluonic correlations; the collinear correlation effects are produced by the tensor 
that multiplies the functions G g a(z) in Eq. (14"T]) . Using Eq. ( I4ip and performing the Fourier 
transformation of Eq. (1571) from q T space to b space, we reproduce the lowest-order structure 
of the resummation formulae in Eqs. ( lH?]) and f l2"U]) (more details are given in Ref. [2H]). The 
term proportional to the unpolarized splitting function P g l a(z) produces the evolution of the 
parton density f g /h(x,fi 2 ) up to the scale /i 2 = b\jb 2 \ this term also produces [12] a residual effect 
(included in the function C ga (z;as) of Eq. fl26]) ) that is related to the definition of the parton 
densities in the MS factorization scheme. The term proportional to Ggl{z) produces the first-order 
contribution to the function G ga (z; as) in Eq. ( 126]) . 



The factorized structure of Eq. ( 120]) has no direct interpretation in terms of a probabilistic 
(classical or quasi-classical) partonic picture. In fact, the structure originates at the cross section 
level from an underlying quantum-mechanical interference between scattering amplitudes and their 
complex conjugates. In particular, the gluonic tensor in Eqs. ( 1201) and ( f26l) (unlike the quark 
coefficient C qa in Eq. ( TT3l) ) cannot be interpreted as the simple residual (and factorization-scheme 
dependent) effect of the customary evolution of the parton density f a /h(x,fi 2 ) (which is a Lorentz 
scalar) from a non-perturbative scale up to the scale /i 2 = b^/b 2 . Incidentally, we note that the 
first-order functions Ggl{z) (a = g, q, q) in Eqs. fl30|) and fl3T|) do not depend on the factorization 
scheme of the parton densities. 

The interference phenomenon that leads to Eq. ( 120]) produces specific physical effects. The 
tensor D^ u (see Eqs. (12T)]) and (|2"g]) ) explicitly depends on the direction of the impact parameter 
vector b in the transverse plane. This dependence produces enhanced (i.e. non-suppressed by 



^This tensor can be replaced by the tensor d^(p,p) of Eq. ([27]) . because of current conservation (gauge invari- 
ance) of the scattering amplitude of the factorized hard-scattering subprocess g{zp) + g{p) — > F. 
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terms of 0(gy/M)) spin and azimuthal correlations in the qx differential cross section da F at 
small values of qx- The spins of the gluons in the gluon fusion process of Eq. f l2~T|) are correlated 
through Eq. f )20|) . The azimuthal angles of the particles of the final-state system F are correlated 
to the azimuthal angle of qx through Eq. ( l20l) and the Fourier transformation in Eq. ( [19]) . 

Note that spin and azimuthal correlations do not necessarily show up their effect simultane- 
ously. For instance, considering the differential cross section averaged over the azimuthal angle 
of qx, the azimuthal correlations cancel, whereas the spin correlations can survive. This obser- 
vation also implies that, in gluon fusion processes, the naive resummation factor of Eq. ffl3|) is 
not correct even in the case of the azimuthally-averaged cross section dopjdq\ (see Sect. [5]). To 
explicitly illustrate this observation in a very simple manner, below we consider a specific process: 
the production of the SM Higgs boson, H, by the gluon fusion mechanism. Since if is a scalar 
particle, its qx cross section has no azimuthal correlations. 

In the case of Higgs boson production, the final-state system F in Eq. (TjQ) is simply H, with 
momentum g M . Within the SM, the corresponding gluon fusion production mechanism is mediated 
by a heavy-quark (mainly, top quark) loop. Our conclusions are unchanged if we consider a generic 
ggH effective coupling (such as, for instance, the SM effective coupling that is obtained in the 
large- M t approximation, M t being the mass of the top quark). Since if is a scalar particle of 
spin 0, the tensor structure of the corresponding hard factor H^~^ mV2 in Eq. (120]) is uniquely 
determined by Lorentz covariance, parity conservation and gauge invariance; we have 



on Mt/M or other parameters of the effective coupling ggH). Using Eqs. (126]) and (PHI) , the Higgs 
boson resummation factor of Eq. (120]) is 



Note that the right-hand side does not depend on the direction of b: this implies that the qx 
distribution has no azimuthal correlations in the small-g^ region. This result is consistent with 
the fact that the q T cross section for H production has no azimuthal correlations at any values 
of qx- 

We can compare the result in Eq. (f45j) with the structure in Eq. ([13]) . Owing to the specific 
factorized dependence on z\ and z 2 of the two terms in the square bracket of Eq. (]4"5]) . the term 
proportional to G(z{) G{z2) cannot be removed by a redefinition of the function C(z;as). This 
shows that the naive gluon fusion expression in Eq. (1131) is not valid even in the simple case of 
Higgs boson production. The structure of Eq. (Hoi) obviously differs from the (correct) structure 
of Eq. ([13]) for qq annihilation processes. The additional term proportional to G{z\) G(z 2 ) is 
indeed produced by (gluon) spin correlations (see also Sect. H]), which have no analogue in qq 
annihilation processes. We finally observe that the term G(z\) G(z 2 ) in Eq. f]45]) starts at 0(a\) 
in QCD perturbation theory (see Eq. f]2"9"]) ). Therefore, its effect first appears in the computation 
of the next-to-next-to-leading order (NNLO) QCD radiative corrections to the Higgs boson q^ 
cross section. 




= HF= H (a s (M 2 )) [C gai ( Zl ;a s (b 2 /b 2 )) C ga2 (z 2 ;a s (b 2 /b 2 )) 
+ Gga^zuasibl/b 2 )) G ga , 2 {z 2 ;a s {bl/b 2 )) } . 



(45) 
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We now comment on our paper in Ref. [IB] . The paper deals with the class of processes in 
Eq. ([1]) and presents a practical formalism for the NNLO QCD calculation of the corresponding 
cross sections at the fully-differential level. The formalism exploits the subtraction method to 
cancel the unphysical infrared (IR) divergences that separately occur in the real and virtual 
radiative corrections. The explicit construction of the subtraction counterterms is based on the 
universal structure of transverse-momentum resummation formulae and on their expansion up 
to NNLO in QCD perturbation theory. These are the resummation formulae discussed in the 
present paper. Working on the research project of Ref. [13], we found the results presented 
in this section, and explicitly documented for the first time in this paper. These results are 
essential for the application of the subtraction method of Ref. [13] to gluon fusion processes. In 
these processes, the naive expression in Eq. (f!3l) does not reproduce the correct (and singular) 
perturbative behaviour of the qr cross section in the limit qx — > 0. Using Eq. f fl3|) . rather than 
the correct result in Eq. (1201) . leads to subtraction counterterms that would spoil the cancellation 
of the IR divergences. Although the theoretical results of the present paper were not explicitly 
illustrated in Ref. [13] . they were actually taken into account in the NNLO computations presented 
therein. In particular, the explicit application to Higgs boson production (which was implemented 
in the Monte Carlo code HNNLO) considered in Ref. [13] is based on and implements the results 
(e.g. Eq. (14"5]) ) illustrated in this section. 

We also add a brief comment on Ref. [28]. The authors of Ref. [28] study transverse- momentum 
cross sections at small by introducing a factorization formalism that differs from the resumma- 
tion formalism considered in Sect. [2] and in this section. In the case of SM Higgs boson production 
or, more generally, gluon fusion processes, a perturbative ingredient of the factorization formulae 
in Refs. [28], [34] is a collinear function (the Lorentz indices fi and v and the parton indices 
g and a refer to the notation used throughout our paper), which emerges from authors' analysis 
based on Soft Collinear Effective Theory. The function I££ is explicitly computed [28| 134"] up to 
its first-order contribution in as- The result at 0(as) is expressed in terms of two form factors, 
J-f a and J-f °, that multiply the tensors in Eqs. (1271) and (]2"8l . respectively. Moreover, at order 
as, we note (see3 Eqs. (35) and (37) in Ref. [34]) that the ratio J-f 9 /J-f 9 is equal to Ca/Cf] the 
same value, Ca/Cf, is obtained by considering the ratio of our functions G gg (z) and Ggq(z) in 
Eqs. (l30j) and (l3~Tj) . In Ref. [28], the authors also note that the form factor J-f does not contribute 
to the Higgs boson qx cross section at the next-to-leading order (NLO). These first-order features 
of the function 1^ have clear analogies with the structure of our universal coefficient tensor 
in Eq. (125]). 



4 The gluon fusion resummation formula in helicity space 

The factorization formula in Eq. (1201 involves sums over the Lorentz indices of the gluons (i.e., of 
the gluon field A^). These sums can be replaced by corresponding sums over the spin polarization 
states of the gluon. In particular, it can be convenient to consider physical polarization states of 
definite helicity A (A = ±). In this section we present the helicity space version of the factorization 
formula in Eq. (1201) . 

tWe assume that there is a typo in the overall sign of the argument of the hypergeometric function o-Fi in 
Eq. (37) of Ref. [34]. 
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Exploiting gauge invariance or, more precisely, current conservation (e.g. px^Cg^^Zi; . . . 
and Eq. (125]) ). the right-hand side of Eq. f )20|) can be written in the following form: 



L x *lgg;a 1 a 2 



Yl i7 (A 1 A),(A 2 ,/ l2 )( x iPi' x 2P2;^;as(M 2 )) 

\±,hi,\ 2 ,h 2 

x C( x a \' h \z i;Pl ,p 2 ,b-,Mbl/b 2 )) C^ 2 M) (z 2 ; Pl ,p 2 M(Xs(b 2 /b 2 )) , (46) 



where A, and hi are helicity space indices (Aj = ±, hi = ±). The relation between the 'helicity 
tensors' in Eq. ( F46|) and the Lorentz tensors in Eq. (120]) is: 



U {\i,hi),(*2,h 2 ) 



(x 1 p 1 ,x 2 p 2 ;^;a s ) = H F ' (xip^x^riias) 



x £ (x l) ( x iPi) l £ (hi)( x iPi)\ £ (x 2 )( x 2P2) [e^ 2 ){x2P2)\ ,(47) 

C^\z l] p 1 ,p 2 ,b;a s ) = [e^\x tPt )Y C^(z i ;p 1 ,p 2 ,b;a s ) e^\x iPi ) , % = 1,2 , (48) 

where e^\p) denotes the polarization vector of a gluon with on-shell momentum p (p 2 = 0) and 
helicity A. 

We recall that the gluon helicity vectors are not uniquely defined. Having chosen the two 
helicity vectors e^\p) and e^\p) = e^\p) , there is still the freedom to change the helicity 
vector basis. For instance, the phase transformation 



e { i\p) e^(p)=e^e^(p) 



i\(p p (A), 



(49) 



defines helicity states eji X \p) that are physically equivalent to e^ X \p) (the phase (p p on the right- 
hand side of Eq. (f4"9l) can also depend on the gluon momentum p) . Note that Eq. ( 1461) is invariant 
under the transformation in Eq. (f4"9l) . whereas the helicity tensors HE h ^ h ^ in Eq. (T47j) and 

in Eq. (1481) are not separately invariant. In any actual computation of HK h ^ , Aa h ^ and 

Cg X a h \ the definition of the helicity basis has to be clearly specified. 



Some general and important properties of the universal coefficient function Cg X d h ^ do not depend 
on the specific definition of the helicity basis. Owing to general relations between helicity vectors 
and using Eq. ( l26l) . we find that the helicity tensor C g d h ^ in Eq. ( 1481) has the following explicit 
structure: 



C^' h ' ) (z t ;p 1 ,p 2 ,b; a s ) = C ga { Zi ;a s ) 6 x » h > + G ga { Zi ; a s ) D {Xl) {p u b) S x ^ 



1,2 , (50) 



where the scalar coefficient functions C ga (z;as) and G ga (z;as) are those of Eq. ( 1261) . 5 x ' h is the 
customary Kronecker symbol (5 + ' + = 5~'~ = 1 , <5 +, ~ = 5~' + = 0), and the helicity coefficients 
D^( Pl ,b), 

2 r . , ^, , n 2 



D^ x \p l ,b) = -- [b-e 



(-A), 



•EiPi 



1,2 



are pure phase factors that have the following explicit form: 



£>W(pi,b) 
D {x) ( P2 ,b) 



3 +2iA [<f>(b)-ipi] 
-2i\ [<P(b)+(p 2 ] 



(51) 



(52) 
(53) 
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Note that the helicity coefficients D^ Xi \pi, b) have a distinctive dependence on 0(b) (which is the 
azimuthal angle of the impact parameter vector b), whereas the phases (pi and y? 2 simply depend 
on the explicit definition (to be specified) of the helicity vectors eft\xipx) and £^(2^2) (see 
Eq. (14~9]) ). As already discussed, the dependence on ipi and <y? 2 cancels in Eq. (H6|) . 

The helicity space representation (150]) of the gluonic coefficient tensors explicitly shows the 
presence of two components: a helicity-conserving component and a helicity-flip component. This 
two-component structure originates from the gluon collinear correlations discussed in Sect. El 
The helicity-conserving component leads to the naive (in the case of gluon fusion processes) 
factorization formula in Eq. (fl3j) . Inserting Eq. (1501) in Eq. (HST) . the helicity-flip component, which 
is proportional to the coefficient function G ga (z; a s ) (see Eqs. (|29|) - (|3T|) ). obviously produces non- 
trivial helicity (spin) correlations in the process-dependent factor H F . 

The helicity-flip coefficients in Eqs. (151Z|) and flBHj) can be rewritten as follows 

£>W(pi,b) = _ e +*M*(b-qT)-w] e +2tt*(q T ) = e + 2iA ^ T ) L>W( Pl ,q T ) , (54) 

£> (A) (p 2 ,b) = _ e -2iA[0(b.q T )+ V2 ] e -2^(q T ) = e -2iA0(b.q T ) £>W q T ) ; (55) 

where 0(qx) is the azimuthal angle of qx (in the centre-of-mass frame of the collision), and 
0(b ■ qx) = 0(b) — 0(qx) is the relative angle between qT and the impact parameter (i.e. b ■ qx = 
bqT cos0(b • qx) )• Using Eqs. (Jl6"j) and fl50l) . the helicity-flip phase factors in Eqs. fl54l and 
(155]) eventually enter the cross section formula in Eq. (1T91) . Since 0(b • qx) is simply the angular 
integration variable of the Fourier transformation in Eq. ( Fl9l) . the dependence on 0(qx) of the qx 
cross section [dap] is directly determined by the phase factors exp (±2iA0(qx)) in Eqs. ( 1541) and 
(|55|) . Therefore, the helicity-flip component of the gluonic helicity tensor in Eq. (|50|) controls both 
the spin correlations and the azimuthal correlations of the qx cross section in the small-gy region. 

In Sect. [3] we have discussed the case of SM Higgs boson production by gluon fusion. The 
corresponding helicity tensor HK^^ , Aa h s is obtained by using Eqs. fj44l) and ff4Tj) . We have: 

^fArJ 1 ),(A 2 ,/ l2 )(^iPi ; ^;«s(M 2 )) = i ^(a s (M 2 )) 5 h ^ e^-^^-Hw) , (56) 

where we have used the relation 

£ ( M A 2 2) (x 2 p 2 ) = -e iAl( ^ 1+ ^ 2) <5 Al ' A2 . (57) 

In Eq. (|56p . the constraint Ai = A 2 (and /ii = ^. 2 ) obviously originates from helicity conservation 
for the production of a boson with spin 0. Inserting Eqs. fl50l . fl52l) . fl53|) and fl56l) in Eq. ff4^1) . we 
reobtain the result in Eq. fj4"5l) . Terms proportional to C(2;i)G(2; 2 ) are absent from Eq. (j35]), since 
helicity conservation in the hard-process factor of Eq. ( l56l) forbids contributions that are produced 
by a single helicity flip. 

We comment on some of the results of Refs. [271 BSJ- Studying diphoton production, the 
authors of Ref. [27] made the following important observation: the description of the small-gr 
behaviour of the diphoton transverse-momentum cross section requires the introduction of new 
logarithmically-enhanced spin-flip contributions (which affect the azimuthal angle dependence of 
the produced diphoton system), through a mechanism that is unique to gluon scattering. This 
general observation was based [27] on the computation of the first-order QCD radiative corrections 



M V Pi ■ P2 J 
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to the lowest-order gluon fusion process gg — > 77. Expanding our resummed formulae (see, e.g., 
Eqs. f|T9|) . fj46|) and f l50|) ) to the first order in as, we find agreement with the structure of the 
diphoton first-order results of Ref. [27] . To be precise, we 'almost' agree with those first-order 
results: our expression contains a term proportional to ^2 a=gqq - G^gl(z) f a / h (x/ z), whereas the 
expression in Ref. [27] includes only the contribution Ggg(z) f g / h (x/z) (the function P g / g {z) of 

Ref. [27] is related to G g l) g (z) in Eq. (ED]); namely, P g/g {z) = ZG g l) g (z) ). The authors of Ref. [27] 
also proposed an all-order generalization of their first-order results. The generalization (which is 
expressed in the Collins-Soper frame [36] of the diphoton pair) is achieved by introducing a new 
helicity-flip unintegrated parton density of the colliding hadrons. This unintegrated (transverse- 
momentum-dependent) parton density, or, more precisely, its Fourier transform to b space is 
denoted by 5 x ~ h V g i h (x, b), where V g , h (x, b) does not depend on the gluon helicity state A = —h. 
The general dependence of V g i h (x, b) on b and the possible relation between V g , h (x, b) and the 
customary parton densities are not specified in Ref. [27J (apart from the relation with f g /h that 
can be inferred from the first-order calculation therein). Our resummation formulae also contain 
one helicity-flip component for each of the colliding hadrons. This component is proportional to 
the following expression: 

f 1 d? 

5 x,- h / _ J- G ga (z;a s (bl/b 2 )) f a/hi (x/z,b 2 /b 2 ) D^( Pi ,b) , * = 1, 2 , (58) 

a=g,q,q 

which is directly related to the parton densities f a /h through the longitudinal-momentum con- 
volution with the perturbatively-calculable gluonic coefficient functions G ga (z;as) (a = g,q,q). 
Note also that the expression in Eq. (|58l) explicitly depends on the gluon helicity A through the 
helicity coefficients D^ x \pi,h) (i = 1,2). The helicity coefficients introduce a definite and non- 
trivial (though, simple) functional dependence on the azimuthal angle 0(b) (see Eqs. (1521 and 
fl53|) ) or, equivalently, on the relative azimuthal angle 0(b • q T ) (see Eqs. fl5ll) and fl55l) ). This 
functional dependence eventually produces definite coherent correlations (see, e.g., Eq. f|45|) and, 
more generally, Eqs. fl66l) . flTUl) and fITBI) in Sects. [5] and [6]) between the helicities, Ai and A 2 , of the 
two colliding gluons. Owing to this helicity dependence of the helicity-flip components in Eq. ( 1581) . 
by using our gluon fusion resummed formulae we are not able to reproduce the diphoton results 
of Ref. [27J (e.g., the structure of Eqs. (27), (30) and (38) in Ref. [27]): the differences show up 
starting from contributions of relative order a| with respect to the lowest-order process gg — > 77. 
Few additional comments about this point are included in the final part of Sect. [61 



5 Azimuthally-averaged cross sections in gluon fusion 
processes 

Starting from the general qT cross section in Eq. (j3J), we can consider its average over 0(qT), at 
fixed values of the additional kinematical variables f2 of the final-state system F. We thus define 
the following azimuthally-averaged cross section: 

/ dap v _ f 27T d0(q T ) da 



d 2 q T dM 2 dydQ'* J 2vr d 2 q T dM 2 dy dfl 

\ dg dM 2 Ty dn (P^«T,M,y,n) . (59) 
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Following the shorthand notation of Eq. OH]), the singular component (in the small-g-r region) of 
the cross section f )59|) is denoted as [dap],. Then we consider the decomposition in Eq. ffl8|) . As 
recalled at the end of Sect. [21 in the case of qq annihilation subprocesses, small-g^ resummation 
does not produce any azimuthal correlations and, therefore, the corresponding [da F ](^ -ann -) anc i 
[dap]^ ann '' ) are obtained by the same resummation formula (see Eqs. ([7|) and ( TlBl with c = q, q). 
In the case of gluon fusion processes, the resummation formula for the azimuthally-averaged cross 
section is obtained from Eq. ( 1T91 by replacing 

(60) 



and by performing the following replacement 



(61) 



in the integrand on the right-hand side; the two-dimensional Fourier transformation con be re- 
placed by the Bessel transformation as in Eq. (fl6|) . Therefore, the resummation formula for the 
azimuthally-averaged cross section is: 



[d(JF 



(g-fua.) 



M 2 



a l5 a 2 21 J** Z2 



d<7 9<l,F 



db 

2^ 



b J {bq T ) S g (M,b) 



f ai /hi(xi/zi,b /b ) fa 2 /h 2 (x 2 /z 2 ,b /b ) , (62) 



where the integrand factor [H F CiC 2 ~\^^ a a2 has the following explicit expression: 

[i^dCj^ = ^(x 1 p 1 ,x 2 p 2 ;f2;as(M 2 ))C 9ai ( 2l ;a s (^/& 2 )) C 9a2 (z 2 ; a s (b 2 /b 2 )) 

+ H F (x 1 p u x 2 p 2 ;n;a s (M 2 ))G gai (z 1 ^s(bl/b 2 )) G ga2 {z 2 ; a s (b 2 /b 2 )) .(63) 

This result can easily be obtained by using Eqs. fj46l) . fl50l . fl54l) and fl55|) and the following 
elementary integral: 



2t 



#(qT) 

2tt 



e ±in</>(q T ) _ fin,0 



71 = 0,1,2, 



(64) 



The process-dependent factor H F on the right-hand side of Eq. (|63|) is defined in Eq. (!24j) . Its 
equivalent representation in helicity space is 



Ai,A2 



The spin-correlated hard-scattering factor if^ has the following expression: 



HG{x 1 pi,x 2 p 2 ;£l;a s ) = z2 H < 



(A-A),(A,-A) 



'x 1 pi,x 2 p 2 ;^l;as) e 



- 2 iX(cpi+lfi2) 



H 



/Ltl I/i ,/12 ^2 



where the 4th-rank tensor d 



■Hi vi,y,2 v 2 
(4) 



IS 



d[^ v \ V ^p 2 ) = - 2 



i ww d 1 ^ 2 + d mi/2 d^ 1 - d^ 2 
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(65) 

(66) 
(67) 

(68) 



with d^ v = d fJ - u (pi,p 2 ) (see Eq. (1271) ). From the helicity space representation ([55]) of the process- 
dependent factor H F , we see that both gluons (with momenta Xipi and X2P2) undergo a helicity 
flip (i.e., Ai = —hi and A 2 = —h 2 ); moreover, we note that the two helicity flips are correlated: 
they occur coherently, with the constraint Ai = A2. 

As already anticipated in Sect. El in gluon fusion processes the naive resummation factor 
of Eq. ( fl3l) is not correct even in the case of azimuthally-averaged cross sections. The correct 
resummation factor in given in Eq. f )63|) . Comparing Eqs. ffl~3]) and ( 163]) . we note the presence 
of the additional gluonic coefficient function G ga (z; as) and of the corresponding hard-scattering 
factor Hq (which, in general, differs from the spin-uncorrelated factor H F ). 



6 Fourier vs. Bessel transformations and azimuthal cor- 
relations of the qx cross section 

Unlike the case of qq annihilation, the gluon fusion resummation factor ^H F CiC 2 \ in Eqs. ( 120]) 
or (I46p depends on the azimuthal angle 0(b) of the impact parameter vector b. Therefore, in 
the resummation formula (I19p we cannot simply use the relation (I16p to perform the azimuthal 
integration involved in the Fourier transformation from b space to q T space. Nonetheless, the 
dependence of [H F C1C2] gg . aia2 011 0(b) is fully specified at arbitrary perturbative orders. More- 
over, this dependence is sufficiently simple to be handled in explicit form. As shown below, the 
azimuthal integration involved in the Fourier transformation can explicitly be performed. The 
two-dimensional Fourier transformation is thus replaced by one-dimensional Bessel transforma- 
tions. The weight functions of these one-dimensional transformations are the Oth-order Bessel 
function Jo(bqr) and higher-order Bessel functions, such as, the 2nd-order and 4th-order functions 
J 2 (bq T ) and J A (bq T ). 

To present the relation between Fourier and Bessel transformations, we first consider the gluon 
fusion resummation formula in helicity space. In this formulation, the 0(b) dependence of the 
integrand factor [H F C 1 C 2 ] (see Eqs. (gg]) and (]50])) is given by the helicity coefficients D^ x \pi, b) 
in Eqs. ([52]) -( 155]) . We have to perform the Fourier transformation of contributions that are linear 
and quadratic with respect to these helicity coefficients. The explicit results of the corresponding 
integration over 0(b) are the following: 

J e ,b ' qT DW( Pil b) F(b 2 ) = DW(p tl q T ) J ^ e ibqT e ±2lX ^ h ^) F(b 2 ) 

r+00 

= -D^\ Pi ,q T ) / dbbJ 2 (bq T ) F(b 2 ) , (69) 
Jo 



•fin the specific case of SM Higgs boson production, we have Hq~ h = Hg~ H (see Eq. (|45j) ) 
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f ^ e ib-q T £,(Al)(p ljb ) D^fab) F(b 2 ) 
J 27T 

= D^\p 1 ,q T )D^\p 2 ,c lT ) [ — e ibqT e +n^)^T) F ^ 

J 27T 
p+oo 

= D^\ Pl ,q T )D^\p 2 ,q T ) / dfe b [5 Al - Aa J (6g T ) + J 4 (bq T ) ] F(b 2 ) 

Jo 

p+oo 

= 5 \ lt \2 e -MAi (vi+va) / db b J {bq T ) F(b 2 ) 

Jo 

p+oo 

+ 5 Al '" A2 D( Al )(pi,qx) ^ (A2) (P2,qT) / dbbJ 4 (bq T )F(b 2 ) , (70) 

where -F(b 2 ) denotes a generic function of b 2 . The results in Eqs. (169]) and (170]) are straightfor- 
wardly obtained by simply using the following integral representation of the Bessel function J 2n : 

J 2 „(x) = (-l) n r p. e +™<x»*e±vn* , n = 0,1,2,... , (71) 

where the variable x is a real number. 

Considering the formulation in terms of Lorentz tensors (see Eqs. ( 120]) and ( 126]) ). the depen- 
dence on 0(b) is produced by the tensor D^ v {pi,p 2 ; b) in Eq. (12 8|). The results analogous to those 
in Eqs. flSSJ and ([70]) are: 

/" r/ 2 h /" +00 
y — e th ^ D^( Pl ,p 2 ;b) F(b 2 ) = - D^( Pl ,p 2 ;q T ) jf dbb J 2 {bq T ) F(b 2 ) , (72) 

y ^ e *-W£,m^(p ljP2;b )£,^(p 1>p2;b ) F(b 2) 
= <r'" 2I >i,P2) / dbbJ (bq T )F(b 2 ) 



(4) 

/•+£» 

+ ^ 4 r i,M2i/2 (pi,P2;qT) / dbbUbq T ) F(b 2 ) , (73) 



(4) 

where d^ Ul,fl2U2 is given in Eq. (168]) . and the 4th-rank tensor n^" 1, ' t2I/a is 

1}™^!^^) = /^^^ . (74) 

The derivation of Eqs. (1721) and (]75j) is similar to the derivation of Eqs. f]6"9~]) and f]70]) . apart from 
few additional algebraic manipulations. 

We note that the three Bessel functions Jo, J2 and J 4 are not independent. Owing to general 
recursion relations between Bessel functions, we have: 

Hx) = 4(6 ~ X2) J 2 (x) - 3 J (x) . (75) 

This relation can be used to express the results in Eqs. ( )6"9l) . fjTU]) . (1721) and (1751) in terms of two 
(rather than three) Bessel functions. 
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The representation in terms of Bessel transformations offers a technical simplification of the 
resummation formula, since the two-dimensional Fourier transformation is replaced by one-dimen- 
sional transformations. Moreover, in Eqs. fl69|) and flTUl) (or, equivalently, in Eqs. f!72|) and f J73|) ) 
the dependence on the azimuthal angle 0(qT?) is fully and explicitly factorized with respect to the 
dependence on the magnitude, qr, of the transverse momentum (the dependence on qp is produced 
by the integration over b). Therefore, we can express the gluon fusion resummation formula in 
a form that manifestly exhibits the functional dependence on 0(qx) of the q T cross section, at 
small values of qx- 



To be explicit, we insert Eqs. (169]) and flTOj) (or, Eqs. (1T21) and 
f[T9~j) . and we obtain: 



3h) in the resummation formula 



[dap 



(g-fus.) 



+ 
+ 



id°F) ClG2 [^ F (0(q T ))] CiG2 + [da F ) GlC2 [# F (<Kq T ))] GiC2 
[da F ) GG [^ F (0(q T ))] ( 



I GG 



(76) 



The notation [if F (0(q T ))] (rather than simply remarks that each of these factors (the 

subscript I is I = C\G 2 ,G\C 2 or GG) depends on q T only through its azimuthal angle 0(qT) 
(i.e., these factors do not depend on the magnitude of qx). All the other terms on the right-hand 
side of Eq. ( 1761) depend on qp, but they are independent of 0(qT)- 



In Eq. ( |76l) . the gluon fusion cross section [dap 



(g-fus.) 



is partitioned into several contributions. 
The first contribution is equal to [dap]^ ius '\ the azimuthally-averaged cross section. Since the 

cross section [Axp]^ -fus '' ) is evaluated at fixed values of 0(qT), our notation is imprecise. The 
notation actually means that the first contribution on the right-hand side of Eq. (!76j) is explicitly 
given by the expression presented on the right-hand side of Eq. (!62l ; this expression coincides 
with the expression of the azimuthal average over 0(qT) of [<iaF]^~ fus '' ) . The other cross section 
contributions, [Atf] 7 , in Eq. (!76j) have the following form: 



[da F ]j 



M 2 
s 



J (°) 



db 

2^ 



b S g {M,b) 



x 



E 



1 dz\ 



1 dz-i 



.t'2 



Z2 



[JCiC^.^ f ai/hl ( Xl / Zl ,bl/b 2 ) f a2/h2 (x 2 /z 2 ,b 2 /b 2 ) , (77) 



where the integrand factors denoted by [JC1C2] 1 are given by the following explicit expressions: 



[JC X C, 
[JC X C, 



1C1G2 

\gg;a 1 a 2 



]GlC2 
\gg\a\ai 

iGG 
\gg;aia 2 



J 2 {bq T ) C gai { Zl -as{bl/b 2 )) G ga2 (z 2 ;a s (b 2 /b 2 )) 

J 2 (bq T ) G gai {zi\as{bl/b 2 )) C ga2 (z 2 ; a s (bl/b 2 )) 
Ji{bq T ) G gai (zi;a s (bl/b 2 )) G ga2 (z 2 ; a s {bl/b 2 )) 



(78) 

(79) 
(80) 
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The hard-scattering factors [iJ F (0(q T ))] , are process dependent; they are 

[# F (0(q T ))] ClG2 = -< 1 ^ 2 , 2 (z 1 p 1 ,z 2 p 2 ;fi;as(M 2 ))d^^(Pi ) P2)^' i2l ' 2 (pi,P2;qT) (81) 

Ai, A2 

= cos(20(q T )) £ iff AiiAiUA2i _ A2) (x lPl ,x 2 p 2 ;fi;«s(M 2 )) e" 2 ^ 2 (82) 
Ai, A2 

- sin(20(q T )) £ zA 2 iff Ai>Ai)i(A2i _ A2) (a; 1 p 1 ,x 2 p 2 ;^a s (M 2 )) e- 2a ^ 2 , 

Ai, A2 

[^ F (0(q T ))] GiC2 = -<, liM2I/2 (x lPl ,x 2 p 2 ;ri;as(M 2 ))D^^(Pi,P2;qT)rf M2l/2 (pi,P2) (83) 
= E <a^a 1 ),(A2,A2)(^iPi^ 2 P 2 ^;«s(M 2 )) e -*M**-««*» 

Ai, A 2 

= cos(20(q T )) £ iJf Ai> _ Ai)i(A2iA2) (x 1 p 1 ,x 2 p 2 ;fi;a s (M 2 )) e" 2iAl ^ (84) 

Ai, A2 

+ sin(20(q T )) ^i^(A 1 ,-A 1 ),(A2,A2)(^i^2P 2 ;^;as(M 2 )) e" 2lAl ^ , 

Ai, A2 

[iJ F (0(q T ))] GG = H^^ 2U2 (x lPl ,x 2 p 2 ;n;a s (M 2 )) D^^( Pl ,p 2 ;q T ) (85) 
= E ^(A,-A),(-A,A)(^i,^ 2 ^;«s(M 2 )) e- 2 ^-^ 2 ^ 

A 

= cos(40(q T )) H( x ^ xu _ x<x) (x lPl ,x 2P2 ;n;a s (M 2 )) e~ 2lX ^~^ (86) 

A 

+ sin(40(q T )) i^H[ Xi _ xu _ XiX) (x lPl ,x 2 p 2] n;a s (M 2 )) e -*M*i-*»> . 

A 

The factors [i/ F (0(qT))j c G and [i? F (^(qT))] GiC2 involve a single helicity flip (see Eqs. ( 1821 
and (1541) ). From the helicity space representation (1561) of the factor [if F (0(q T ))] , we see that 
both gluons (with momenta X\p± and x 2 p 2 ) undergo a helicity flip, and the two helicity flips are 
correlatedby the constraint Ai = — A 2 . We recall that the double helicity flip with Ai = A 2 leads to 
the factor Hq (see Eq. (166]) ) that enters the cross section contribution [<icr i r]^~ fus '- ) (see Eq. (163]) ). 



By direct inspection of Eqs. (176]) . (152"]) . (154]) and (1561) . we note that the 0(qi) azimuthal de- 
pendence of the logarithmically-enhanced terms at small qp is fully determined by the general 
structure of the transverse- momentum resummation formula. The qT cross section [dap ]( 9-fus -) 
contains a contribution that is independent of 0(q-r) (namely, the term [dap]^ fus '' ) ) plus a lin- 
ear combination of the four trigonometric functions cos(20(q T )), sin(20(q T )), cos(40(q T )) and 
sin (40(q T )) . No other functional dependence on 0(qx) is allowed by the gluon fusion resummation 
formula. 

A general remark about the azimuthal dependence is required. Since we are dealing with colli- 
sions of spin unpolarized hadrons, the corresponding cross sections are invariant under azimuthal 
rotations in the transverse plane. Therefore, the multidifferential cross section [<icrp] ( ' 9 ~ fus '' ) cannot 
depend on the absolute value of (fiiq?)- The azimuthal dependence of [dap]^ 9 ~ ius '^ can only appear 
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through final-state azimuthal correlations, namely, through functions of relative azimuthal angles 
A0j, such as, for instance, A0j = 0(q*r) — <KqTj)- Here, 0(qT«) denotes the azimuthal angle 
of one of the particles in the produced final-state system F (see Eq. ([]])). Note that these az- 
imuthal correlations are consistent with our previous conclusions about the functional dependence 
on 0(q-r); indeed, we have cos(2A0j) = cos(20(q T )) cos(20(q Ti )) + sin(20(q T )) sin(20(q Ti )) , 
and analogous relations apply to sin(2A0j), cos(4A0j) and sin(4A0j). According to the general 
notation that we have used in this paper, the 0(qTi) dependence of the multidifferential cross 
section in Eq. (j3J) (or, Eq. (jSJ)) is introduced through the final-state kinematical variables gener- 
ically denoted by f2 = {0^, Ob, • • • }. In Eqs. (1521) . (!5i|) and (!86|) . the dependence on Q (e.g., 
on 0(q Ti )) of H F (xiPi,X2P2', Q Ois{M 2 )) combines with the explicit dependence on 0(q-r) to pro- 
duce the final-state azimuthal correlations. The functional form of the azimuthal correlations is 
determined by the 0(qT?) dependence of the hard-scattering factors [if F (0(qx))] v through per- 
turbative coefficients (we recall that the Lorentz tensor H F Ui M2 U2 , and the equivalent helicity 
tensor HK h ^ , Aa h ,, are power series functions of as) that are process dependent and observable 
dependent (i.e., they depend on the physical observable that is specified by the multidifferential 
cross section dap/dfl). In particular, if the multidifferential cross section is insensitive to the 
azimuthal angles 0(qTi) (see, e.g., the simple case of inclusive production of the SM Higgs boson), 
the corresponding hard-scattering factors [if F (0(q T ))] 7 vanish order-by-order in perturbation 
theory. 

The structure of the gluon fusion resummation formula (1761) is much richer than the structure 
of the corresponding 'naive' (i.e., extrapolated from qq annihilation) formula in Eq. ([7j). The ad- 
ditional structure is due to the helicity-flip contributions. These contributions are perturbatively 
driven by the gluonic coefficient function G ga (z; as), and they lead to several hard-scattering 
factors: the factor H F (see Eq. (166]) ). which contributes to the term [daF^f f , and the fac- 
tors [if F (0(q T ))] . We briefly sketch the small-g^ singular behaviour produced by these various 
terms in QCD perturbation theory. To this purpose, we expand the integrand of the resummation 
formula in powers of as = as(n%) (with ji R ~ M), and we explicitly perform the Bessel transfor- 
mations from b space to qx space (technical details on this procedure, and explicit perturbative 
formulae can be found, for instance, in Ref. 



The single helicity-flip terms first contribute at the relative order as- At this perturbative 
order, they lead to a partonic cross section contribution that is proportional to 



O's 



7T 



5{l- Zl )5 gai Gfl 2 {z 2 ) [H F (<j>( qT ))} 



C1G2 



57) 



and to an analogous contribution that is obtained by the exchange 1 -H- 2 of the subscripts. Since 
at this order there is a leading logarithmic term of the type 4rhi(M 2 / g|,) (this term is due 

\_1t J + 

to the 'naive' contributions, with no helicity flips), the term in the expression (!57|) represents a 
next-to-leading logarithmic effect. The double helicity-flip term first contributes at the relative 
order a§; it produces a partonic cross section contribution that is proportional to 



as 

7T 



GG 



where the functions Ggl{z) are given in Eqs. (15UI) and (15T1) . 



q\ 
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As we have discussed in Sect. HI our gluon fusion resummation formula produces differences 
with respect to the diphoton results presented in Ref. [27]. These differences are evident starting 
from contributions at the relative order a§. To explicitly point out the 0(a\) differences, we can 
rewrite the structure of Eq. (38) in Ref. [27] by using our notation: this gives an expression that 
is proportional to 



osy 

7T / 



G<fi{ Zl ) G g % 2 ) 5 gai 5 ga , 2 2 ( [# F (0(q T ))] 



I GG 





" 1 " 


) 


.Qt. 



39) 



where the sum of hard-scattering factors in the round bracket originates from Eq. (30) of Ref. 
The qx dependence of the expression (18"9"j) differs from that of our expression (jSSl) . and the difference 
is not removed by setting Ggq(z) = G g q(z) = (the difference between the 0(as) expressions in 
Eq. (33) of Ref. [27] and in our Eq. fl57|) disappears by forcing G^^z) and G g q(z) to vanish). 



At higher perturbative orders, the small-g^ singular behaviour in the expressions (jSTJ and 
is enhanced by powers of ln(M 2 / q 2 ^) . The dominant logarithmic enhancement is produced by the 



gluon form factor S g (M,b), which appear in [dap 



(s-fus.) 



'see Eq. (1621) ) and in each cross section 



contribution [d0>]j (see Eq. (1771) ). As is well known, at the relative order ctg, the customary (i.e., 
with no helicity-flip contributions) leading logarithmic terms have the following structure: 



I In— 
L Qt \Qt 



n > 1 



(90) 



At this order, the gluon fusion helicity-flip contributions produce the following logarithmic be 
haviour: 

1 hM f M 2 



[# F (<Kq T ))] 



C 1 G 2 



GG 



In 

(It 



2n-5 



4" In 

Qt 



M 2 
Qt 



Qt 



+ 



2n-4 



+ 



M 2 
Qt 



+ 



n > 2 



n > 3 , 



n > 3 . 



(91) 
(92) 
(93) 



The dots in the round brackets of Eqs. ( T90T) (T93T) stand for subdominant terms in each corre- 
sponding expression. The comparison between the expressions fl90|) and fl9TT) shows that the single 
helicity-flip terms produce next-to-leading logarithmic contributions at each perturbative order. 
The double helicity-flip terms lead to subdominant logarithmic contributions. 



The singular behaviour that is observed order-by-order in the QCD perturbative expansion 
is cured by qT resummation. We recall that the resummed gluon form factor S g (M,b) (and, 
analogously, the quark form factor S q (M, b) in the qq annihilation channel) provides the integration 
over b in Eqs. fl62l) and flTTj) with a strong damping factor in the large-6 region (roughly speaking, in 
the region where b > 0(1/ M)). This damping effect (the simple resummation of the leading double- 
logarithmic terms, \n 2n (b 2 M 2 ), in b space is sufficient to highlight the effect [2] ) eventually leads 
to resummed perturbative predictions for the qx cross section that are physically well-behaved 
in the small-g^ region. In particular, the qualitative behaviour of the resummed qx cross section 
(176]) at very low values of qT can be examined by performing the limit qT — > of Eqs. (|62l) and 
(J77J. In this limit, we can write: 



[dap 



(g-fus.) 



[dap 



(qq—ann.) 



const. 



(94) 
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l da F} Cl G 2 ~ l d(T F} Gl c 2 ~ ?T > ( 95 ) 

[oM GG ~ g * . (96) 

The constant behaviour in Eq. (194"]) is just a result of Ref. [2j. It follows [2] from a simple reasoning 
(and a minor modelling of the Sudakov form factor at very large values of b, b ~ 0(1/ Aqcd) )• In 
few words, the reasoning amounts to the observation that the presence of the resummed Sudakov 
form factor justifies the use of the low-g^ approximation Jo(^t) ~ 1 to extract the behaviour of 
the resummation formula fl62l) at qx ~ 0. The behaviour in Eqs. fl95|) and ( 196]) follows from the 
same reasoning. We simply note that, in the case of the helicity-flip components [dap]j of the 
qx cross section, the resummation formula ( 1771) involves higher-order Bessel functions. Thus, we 
have just used the \ow-qx approximation J2(bqr) ~ b 2 q£ for the single helicity-flip components 
(see Eqs. ( 178]) and ( 179]) ) and the corresponding approximation J^bqx) ~ b A q£ for the double 
helicity-flip component (see Eqs. (180]) ). 



7 Summary 

Considering the hard-scattering production of high-mass systems in hadron-hadron collisions, in 
this paper we have examined the corresponding transverse-momentum cross sections at small val- 
ues of q T . We have presented a study of the contributions that are logarithmically enhanced 
order-by-order in QCD perturbation theory. The enhanced contributions have the form of sin- 
gular (^-distributions of the type 4- ln n (M 2 /gf,) . The all-order analysis and the perturbative 

\_1t J + 

resummation of these terms, in processes (such as the DY process) that are controlled by the qq 
annihilation channel, is a classical QCD result. We have pointed out that this result does not 
suffice for the treatment of processes (such as SM Higgs boson production) that are controlled by 
(or, simply, receive contributions from) the gluon fusion channel. The difference between the qq 
annihilation and gluon fusion channels originates from correlations that are intrinsically related 
to the collinear 'evolution' of the colliding hadrons into gluon partonic states. 

We briefly summarize the main features of our general results on q? resummation in gluon 
fusion processes. The gluon fusion resummation formula for generic qx cross sections is pre- 
sented in Eqs. ([19]) . (120]) and (126]) . The resummation formula controls all the singular (and 
logarithmically-enhanced) perturbative contributions to the qx cross section in the small-g^ re- 
gion. Gluon collinear correlations produce new (with respect to qq annihilation) structures from 
the perturbative evolution of the parton densities of the colliding hadrons. The additional struc- 
ture (see Eq. (126]) ) enters qr resummation through the factorization formula (1201 . The terms 
due to collinear correlations lead, in general, to next-to-leading logarithmic contributions to the 
qx cross section (the leading logarithmic contributions still come from soft-radiation effects in- 
cluded in the customary Sudakov form factor). Gluon collinear correlations are directly related 
to helicity-flip phenomena in the gluon fusion hard-scattering subprocess (see Eq. (I50p ). These 
spin correlations originate from a quantum-mechanical interference: the flip occurs between the 
gluon helicity states in the scattering amplitude and in the complex-conjugate amplitude. The 
helicity-flip phenomenon due to gluon collinear correlations leads to definite correlations between 
the azimuthal angles of the particles in the high-mass system that is produced by the gluon fusion 
mechanism. These azimuthal-correlation effects accompany the dominant (singular) q? behaviour 
of the perturbative cross section in the small-gr region (azimuthal correlations with a similar qx 
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behaviour are absent if the high-mass system is produced by qq annihilation). The functional form 
of the azimuthal correlations is fully specified, at arbitrary perturbative orders, by the gluon fusion 
resummation formula (see Eq. f!76|) and Eqs. ( IHTji -f lHB"]) ). The double helicity-flip component of the 
qx cross section is characterized by a coherent interference between the spin-flipping gluons from 
the two colliding hadrons (see Eq. (1701)). As a consequence of this interference, the double helicity- 
flip contribution produces two distinct terms with a different gr behavior: a term with azimuthal 
correlations (the last term on the right-hand side of Eq. (!76l) ) and a term, with no azimuthal 
correlations, that also contributes to azimuthally-averaged cross sections (see Eqs. (|62|) and (163|) ). 
Gluon collinear correlations thus imply that the differences in the structure of gr resummation 
between the gluon fusion and qq annihilation channels persist even after having performed the 
integration over the azimuthal angle of the transverse-momentum vector. 

An interesting and relevant issue regards the extension of qx resummation to processes whose 
final-state system F contains strongly-interacting particles (partons) such as, for instance, high- 
p± hadrons and, more generally, jets. The general extension to this type of final-state systems 
(which have not been considered in this paper) is still lacking. It requires a proper treatment of 
soft radiation [37] in multiparton hard scattering, namely, the hard scattering of the two colliding 
partons and the final-state QCD partons in the system F. The main features of collinear radiation 
from the two colliding partons are not affected by the presence of the additional hard partons in 
the final state. Therefore the structure of the gluon collinear correlations that we have found 
and documented in this paper is relevant to any extensions of qx resummation. In particular, the 
convolution of the parton densities with the perturbative gluonic tensor in Eqs. f[2"6l) or floUl) is 
expected to appear in the qx resummation formulae for the production of final-state systems that 
contain colour-charged partons. Note that these systems can be produced by qq annihilation (or, 
generally, qq and qq scattering), gluon fusion (or, generally, gg scattering) and gq(q) scattering 
subprocesses on equal footing. Therefore, there is no escape from (almost) ubiquitous collinear 
correlations due to initial-state gluon hard scattering. 
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